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ABSTRACT 

We present Cosmic Microwave Background (CMB) maps from the Santa Barbara HACME balloon 
experiment (Staren et at 2000), covering about 1150 square degrees split between two regions in the 
northern sky, near the stars 7 Ursae Minoris and a Leonis, respectively. The FWHM of the beam 
is ~ 0.77° in three frequency bands centered on 39, 41 and 43 GHz. The results demonstrate that 
the thoroughly interconnected scan strategy employed allows efficient removal of l//-noise and slightly 
variable scan-synchronous offsets. The maps display no striping, and the noise correlations are found 
to be virtually isotropic, decaying on an angular scale ^ 1°. The noise performance of the experiment 
resulted in an upper limit on CMB anisotropy. However, our results demonstrate that atmospheric 
contamination and other systematics resulting from the circular scanning strategy can be accurately 
controlled, and bodes well for the planned follow-up experiments BEAST and ACE, since they show 
that even with the overly cautious assumption that l//-noise and offsets will be as dominant as for 
HACME, the problems they pose can be readily overcome with the mapmaking algorithm discussed. 
Our prewhitened notch-filter algorithm for destriping and offset removal should be useful also for other 
balloon- and ground-based experiments whose scan strategies involve substantial interleaving. 



1. INTRODUCTION 

One of the main goals of the next generation of Cosmic 
Microwave Background (CMB) experiments is to produce 
maps that can resolve degree scale features correspond- 
ing to acoustic peaks in the angular power spectrum (see 
e.g. Bond 1996 and Hu et al. 1997 for reviews), since this 
potentially allows accurate determination of cosmological 
parameters (Jungman et al. 1996; Bond et al. 1997; Zal- 
darriaga et al. 1997). 

To date, maps with degree scale angular resolution 
have been pubHshed from the ACME/MAX (White & 
Bunn 1995), Saskatoon (Tegmark et al. 1997 - hereafter 
"T97d"), CAT (Jones 1997) and QMAP (Devlin et al. 
1998; Herbig et al. 1998; de Oliveira-Costa et al. 1998) 
and Python (Coble et al. 1999) experiments. The maps we 
present are based on the 1996 flight of the HACME bal- 
loon experiment, which is described in detail by Staren et 
al. (2000, hereafter "S2000"). They cover regions around 
the stars 7 Ursae Minoris (7UMi) and a Leonis (aLeo), as 
shown in Figure |l| (the COBE DMR map is from Bennett 
et al. 1996). The data analysis we describe below is an 
interesting pathfinder for upcoming mapping experiments 
since HACME differed from its predecessors ACME/MAX, 
Saskatoon and CAT in three important ways: 

1. The time-ordered data set is substantially larger 
(consisting of ~ 10 temperature measurements). 



2. l//-noise plays a significant role. 

3. The scan strategy interconnects the observed pixels 
in a complicated way. 

To use this interconnectedness for efficient l//-noise 

removal, we apply the lossless mapmaking algorithm de- 
scribed in Tegmark (1997a, hereafter "T97a") tailored to 
the noise power spectrum at hand. Our analysis provides 
the first large scale test of this method, since the large size 
of our data set forces us to use all the numerical tricks 
that have been proposed for accelerating it (Wright 1996; 
Tegmark 1997c - hereafter "T97c"). We describe our anal- 
ysis in Section 2 and present the results in Section 3. 

2. METHOD 

We make our maps using the minimum- variance method 
described by Wright (1996), T97a and T97c, which we 
summarize briefiy below. 

Let Ui denote the temperature measured by the HACME 
experiment at the i*^ observation (there were 312.5 such 
measurements per second), and group this time-ordered 
data set (TOD) into an M-dimensional vector y. As de- 
scribed in S2000, there are M = 3, 323, 250 measurements 
in the 7UMi region and 1,743,625 in the aLeo region, re- 
spectively, for each of the three frequency channels. Dis- 
cretizing the map into N pixels centered at unit vectors 
fi, ...,f AT, we let Xi denote the true beam-smoothed CMB 
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Fig. 1. — Location of observing regions in galactic coordinates. 



temperature in the direction f ^ and group these tempera- 
tures into an A'^-dimensional map vector x. The TOD is 
related to the map by 

y = Ax + n, (1) 

where n denotes a random noise vector and the M x N 
matrix A encodes the HACME scan strategy. Letting Ni 
denote the number of the pixel pointed to at the i*-^ ob- 
servation, we have Ay = 1 if TV^ = j, Aij — otherwise, 
so all the entries of A are zero except that there is a sin- 
gle "1" in each row. We find that n can be accurately 
modeled as a multivariate Gaussian random variable with 
zero mean ((n) =0), so it is completely specified by its 
covariance matrix N = (nn*). 

The problem at hand is to compute a map x that es- 
timates the true map x using the TOD y. As shown in 
T97a, all the cosmological information from the TOD is 
retained in the map x given by 

i= [A*MA]-iA*My, (2) 

where the matrix M = N^^. Substituting this into equa- 
tion (|l|) shows that we can write x = x + e, i.e., that our 
map estimate x equals the true map x plus a Gaussian 
noise term e that is independent of x. This holds for any 
choice of M. This pixel noise has zero mean ((e) = 0) and 
a covariance matrix given by 

£ = (££*} = [A*MA]"i[A*N"iA][A*MA]"\ (3) 

To minimize the noise, we want to choose M k, (at 
least approximately), which gives S near the best possi- 
ble case [A*N~^A]~^. For the simple case of white noise 
from the detector, we would have N cx I, the identity ma- 
trix, and X would be obtained by simply averaging the 



measurements of each pixel in the map. We will refer to 
this as the straight averaging mapmaking method, in con- 
trast to the prescription of equation (|^) , which we will call 
the minimum noise method. In our case, the matrix N 
is far from diagonal, since long-term 1// drifts introduce 
correlations between the noise Ui at different times. To 
eliminate the need for a numerically unfeasible inversion 
of N in equations (H) and , we employ the prewhitening 
trick of Wright (1996) and define a new data set y = Dy, 
where the prewhitening matrix D is chosen so that the 
prewhitened noise n = Dn becomes as close to possible 
to white noise, i.e., so that the filtered noise covariance 
matrix N = (nii*) = DND* « I. Defining A = DA, 
this allows us to rewrite equation (^ as y = Ax + n, so 
equations (^ and (||) give 

i= [A* MA] -1 A* My 

= [A*D*MDA]-iA*D*MDy, (4) 
S = [A*D'MDA]-i[A*D*MNM*DA][A*D*MDA]-\ 

where we make the choice M = I since N w I. To ac- 
celerate the numerical evaluation of these expressions, we 
take advantage of the fact that D*D is almost a circu- 
lant matrix, as described in detail in T97c. Since the 
statistical properties of the noise are approximately con- 
stant over time (the covariance (niUj) depends only on 
the time interval between the two observations), the 

prewhitening filter D is simply a convolution filter. We 
compute it as described in T97a, and the result is plot- 
ted in Figure ^ The prewhitening procedure renders the 
data insensitive to the mean (monopole) of the map. This 
makes the matrix [A*MA] non-invertible, which is reme- 
died using the pseudo-inverse technique described in the 
Appendix of Tegmark (1997b). 

In addition to the map x, we also compute the Wiener 
filtered map x^, which is given by (Wiener 1949) 

x^ = S[S-hN]-ix, (5) 

where the CMB signal covariance matrix S = (xx*) is 
given by the standard formula 

S^,^f2^^^Pdh-r,)B^C^ (6) 

i=2 

for some assumed fiducial power spectrum Ci. In our 
analysis, we approximate the HACME beam as a circu- 
lar FWHM= 0.77° Gaussian, although the actual beam 
has an asymmetry of order 6%. Thus Bi w eit+i)/2^ 
where 6 = (81n2)-i/2 FWHM. As discussed in T97a, the 
Wiener filtered map x^, contains the same cosmological in- 
formation as X (which we will refer to as the unsmoothed 
map to distinguish it from x^,), but is more useful for vi- 
sual inspection since it is less noisy. Wiener filtering has 
been successfully applied to both CMB maps {e.g., Bunn 
et al. 1994, 1996; T97d) and galaxy surveys {e.g.. Fisher 
et al. 1995; Lahav et al. 1994; Zaroubi et al. 1995). 

After removal of calibration data from z, segments dur- 
ing which control commands were transmitted to the gon- 
dola are omitted, which reduces the size of the TOD by 
about 0.4%. The remainder has scan-synchronous offsets 





Fig. 4. — Wiener-filtered maps of the 7UMi (top) and aLeo (bottom) regions at 39 GHz (left), 41 GHz (center) and 43 GHz (right). 
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described in S2000. We determine the time-independent 
scan- synchronous component simultaneously with the 
map by including unknown offsets, one corresponding to 
each of the 125 sampling positions along the circular scan, 
as 125 additional "pixels" in the vector x and widening A 
accordingly with an additional "1" in each row. This trick 
was subsequently applied also in the analysis of QMAP — 
see de Oliveira-Costa et al. (1998). 
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Fig. 2. — The prewhitening procedure is illustrated in the fre- 
quency domain (left) and in the time domain (right) for the 7UMi 
43 GHz channel. The noise power spectrum (upper left) contains a 
l//-component which causes the noise to be almost perfectly cor- 
related between measurements close together in time (upper right). 
By multiplying the Fourier transformed signal by an appropriate 
prewhitening filter (middle left), which corresponds to applying a 
convolution filter (middle right) to the signal, the filtered TOD y 
obtains a noise power spectrum which is close to white noise with 
notches at Hz and the spin frequency (lower left), corresponding 
to a time autocorrelation function that falls of rapidly (lower right). 




Fig. 5. — Noise properties of the unsmoothed (top) and Wiener 
filtered (bottom) 7UMi maps at 43 GHz. The leftmost maps show 
the rms noise in each pixel. The remaining maps show the corre- 
lations with two typical pixels, located at the centers of the white 
dots. 



3. RESULTS 

3.1. 1/f-noise removal 

Figure || (left panel) shows the 43 GHz aLeo map ob- 
tained with the straight averaging method (simply aver- 
aging the observations of each pixel), and illustrates why 
the more sophisticated minimum noise method (middle 
panel) is needed. The visible striping along the circu- 
lar scans (see S2000 for details about the scan pattern) 
is caused by l//-noise, which creates strong correlations 
between map pixels that are scanned within a short time 
of one another. The middle panel shows that our improved 
mapmaking technique has eliminated this striping. This is 
more vividly illustrated in Figure |^ (upper right panels) , 
showing two typical rows of the noise covariance matrix 
S. Since row i gives the correlations to i*'' pixel, the 
small round white dots tell us that strong noise correla- 
tions persist only between pixels that are close neighbors. 
Moreover, these correlations are seen to be quite isotropic, 
with no signs of elongation along scan directions (the only 
exceptions are pixels close to edges). There is a slight 
anticorrelation to very distant pixels caused by the above- 
mentioned monopole removal, which forces the correlation 
map to have a vanishing average. This removal of 1/ f- noise 
also causes a substantial reduction in the rms pixel noise 
(upper left panel in Figure 5), which is about a factor of 
two lower than for the straight averaging map. Both of 
these desirable properties of the pixel noise can be read- 
ily understood by comparing our scan strategy to the four 
toy models described in T97c. Since the circular HACME 
scans cross one another at virtually every pixel, the ge- 
ometry is similar to the "fence pattern" , which was found 
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to produce isotropic noise correlations like those in Fig- 
ure 1^. A "fence" map made with the straight averaging 
method would correlate the noise in a pixel with neighbors 
along the two perpendicular scan directions, which form a 
"+" -shaped region around the pixel. With the minimum- 
noise method, on the other hand, a pixel is linked to all its 
neighbors via a sequence of intersecting scans in a zig-zag 
pattern. This is accomplished automatically in the matrix 
inversion step, and has the effect of isotropizing correla- 
tions (thereby reducing long-range correlations) and utiliz- 
ing the available information better (reducing the resulting 
pixel noise). 

3.2. Wiener-filtered maps 



The rms pixel noise <Ti 



,1/2 



in the unsmoothed 43 



GHz aLeo map (middle panel of Figure ^ is about 400 fiK 
per (0.3°)^ pixel in the central parts and 250 /iK when re- 
binned into (0.8°)^ pixels, i.e., substantially larger than 
any expected CMB signal. What the eye perceives is 
therefore mostly noise, so it is desirable to increase the 
signal-to-noise ratio by smoothing. Moreover, some pixels 
near the edges are extremely noisy (with an rms around 
10 mK) since they were only observed a small number of 
times (upper left panel of Figure |^). To prevent them 
from dominating the image, it is desirable to give them 
less weight. The Wiener filtered map (right panel of Fig- 
ure H) automatically remedies both of these problems. The 
new noise map (lower left panel of Figure g) shows that 
the noise explosion near the edges has been eliminated, 
and the new correlation maps (lower right panels of Fig- 
ure H) show that smoothing has widened the correlated 
region. To avoid imprinting features on any particular 
angular scale, we have used a flat fiducial power spectrum 
Cf = 247rQV5^(^ + 1) in equation (|) just as in T97d, nor- 
malized to Q = 30/xK. Figure |J shows the Wiener-filtered 
7UMi and aLeo maps for all three channels, on a sym- 
metric linear grey scale from low (black) to high (white) 
temperatures. Since some maps are substantially noisier 
than others and the Wiener filtering suppresses their tem- 
perature scale accordingly, the grey scale has been sepa- 
rately normalized for each map to avoid saturation and 
allow direct comparison of spatial morphology. 

3.3. Verification and testing 

We tested our data-analysis pipeline by generating mock 
sky maps that were observed with the actual HACME scan 
pattern. These simulated TOD sets were then run through 
the pipeline, recovering the original maps. When adding 
noise to the TOD with the measured HACME power spec- 
trum, the pixel noise in the recovered map was found to be 
fully consistent with the noise covariance matrix S com- 
puted by the pipeline. 

To test if the reconstructed map was sensitive to the 
pixel size, the analysis was repeated for a range of pixel 
sizes. As expected, this produced virtually identical maps, 
as long as the pixel size was smaller than the 0.77° beam 
size by at least the Shannon oversampling factor 2.5. 

To test if the Wiener-filtered maps were sensitive to the 
choice of fiducial band power, they were computed for the 
band power normalizations Q = 20, 30 and 40 ^K. The 
visual difference between the two extreme cases was min- 
imal, corresponding merely to an increase in the effective 



smoothing scale by about 10% (for more details on why 
this is expected, see the discussion in T97d). 

3.4. Constraints on systematics and the CMB power 
spectrum 

The initial version of our analysis assumed that the scan- 
synchronous offsets were time- independent. The result 
was a marginal detection of fiuctuations in the most sensi- 
tive channel (43 GHz) and upper limits at 39 and 41 GHz. 
However, subsequent tests showed that variations in these 
offsets at the level of a few percent could masquerade as 
signal at this level. Since the offsets are completely dom- 
inated by the first harmonic (at the scan frequency), we 
repeated the entire analysis after throwing away all the 
information about this frequency. This corresponded to 
replacing the whitening filter D discussed in Section |^ by 
a notch filter that vanished at this frequency (as well as at 
Hz). The Wiener- filtered maps shown in Figures I and I 
were made without the notch filter, whereas those in Fig- 
ure I used it. The result of the notch filter was an increase 
in the rms noise of about 50% for the least noisy pixels in 
the raw maps, as well as a reduction of the signal seen 
to levels that were no longer statistically significant. The 
pixel noise ct^ (from detector noise only) is 385 fiK in the 
best pixel in the 43 GHz map of the 7UMi region. The 41 
and 39 GHz channels are noisier by factors of 1.2 and 1.9, 
respectively. For comparison, a CDM-like power spectrum 
normalized to a band power Qflat = 30, fj,K would give 
corresponding pixel fluctuations of order 80 fiK in the this 
map including the mean removal. Most of the fluctuations 
seen in the raw maps are therefore noise rather than sky 
signal, and as described below, this remains true for the 
Wiener-filtered maps as well. The reason that many fea- 
tures are seen to recur between the channels is that their 
noise is correlated. Although there is a slight sensitivity 
to the dipole (with sky rotation marginally breaking the 
degeneracy between the scan-synchronous offset and the 
dipole), the resulting noise error bars on this mode are 
too large to allow a detection after our notch filtering. 

To place quantitative constraints on the systematics and 
sky signal in the maps, we perform a signal-to-noise eigen- 
mode analysis (Bond 1994; Bunn & Sugiyama 1995). We 
define a new data vector 



B'x, 



(7) 



where b, the columns of the matrix B, are the eigenvectors 
of the generalized eigenvalue problem 



Sb = ASb, 



(8) 



sorted from highest to lowest eigenvalue A and normalized 
so that b*Nb = I. A sample of eigenmodes b are shown in 
Figure |^, and the corresponding expansion coefficients yi 
are plotted in Figure ^ If there are no systematic errors or 
residual offset problems, detector noise should contribute 
uncorrelated fiuctuations of unit variance to each mode. 
All systematic problems we have considered would add 
rather than subtract map power. Figure ^ indeed shows 
no evidence of outliers or other problems with our noise 
model. 

The covariance of the expansion coefficients is {yiyj) = 
Sij{l -|- Ai), so cosmic signal at the assumed level would 
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boost the variance somewhat. The very best 43 GHz 7UMi 
mode {i — 1) has A; « 0.32, corresponding to an expected 
rms of {l + XiY^^ ~ foi' Vi- Figure|7|is seen to be quite 
consistent with this. Thus no single mode has signal-to- 
noise exceeding unity. 



100 aOO 300 400 500 

S/N mode number i 

Fig. 7. — The triangles show the first 500 S/N eigcnmode ex- 
pansion coefficients yi for the 43 GHz 7UMi map. If there were 
no systematics or CMB fluctuations in the map, merely detector 
noise, they would have unit variance, and about 68% of them would 
would be expected to lie within the yellow/light grey band. For 
a ffat Qflat = 30 CMB power spectrum, the standard devia- 
tion would be larger for the best modes, as indicated by the thin 
blue/grey band on the left side. The wiggly green/grey curve is the 
rms of the data points yi, averaged in bands of width 25. 

To place cosmological constraints, we need to com- 
bine the information from all modes. We do this using 
the decorrelated quadratic estimator method described 
in Tegmark (1997b). We find that no channel gives a 
statistically significant detection. The sharpest upper 
limit comes from the 43 GHz 7UMi data, giving ST = 
[£{£ + l)C£/27r]i/2 < 77 at 95% confidence with the 
window function shown in Figure 0. This corresponds to 
the effective angular scale i = 38^^. 

What was the price paid for using the notch filter? Fig- 
ure || shows that this destroys information on scales ^ < 20, 
so the 50% increase in rms pixel noise that we found above 
is caused largely by increased detector noise fluctuations 
on these large angular scales, comparable to the size of the 
basic 10° scanning circle. 

4. CONCLUSIONS 

We have presented CMB maps from the HACME bal- 
loon experiment, which was devised as a pathfinder for 
the long-duration balloon experiments BEAST and ACE. 
Our analysis has clearly demonstrated that neither sub- 
stantial l//-noise nor slightly variable scan-synchronous 
offsets are "show-stoppers", but can be efficiently dealt 
with using the matrix-based minimum variance notch fil- 
ter technique. However, to provide interesting cosmolog- 
ical constraints, the raw sensitivity needs to be substan- 
tially improved over the HACME pathfinder specifications. 
Comparing this with the straight averaging method (av- 
eraging the observations of each pixel) , we found that our 
analysis technique improved the results in two ways: 




Fig. 8. — The window functions are shown for the analysis with 
(solid) and without (dashed) the notch filter. The latter is seen to 
destroy sensitivity on scales £ < 20. 



1. The resulting rms noise in the maps is reduced by 
about a factor of two, the equivalent of four times 
more flight time. 

2. Long-range anisotropic noise correlations ("strip- 
ing") is virtually eliminated from the maps. 

This bodes well for upcoming spin-chopped total power 
experiments such as BEAST and ACE. BEAST is antici- 
pated to have nearly 3 times higher sensitivity in fiK s^/^ 
and 2-6 times the angular resolution, with the same ratio 
of 1// knee to spin rate as HACME and 50 times more in- 
tegration time. Indeed, the notch filter approach should be 
useful for any ground- or balloon-based CMB experiment 
using an interconnected scan-strategy. 

Support for this work was provided by NASA through 
grants NAG5-6034 and NAGW-1062 and through a Hub- 
ble Fellowship, HF-01084.01-96A, awarded by the Space 
Telescope Science Institute, which is operated by AURA, 
Inc. under NASA contract NAS5-26555, and by the NSF 
Center for Particle Astrophysics. JWS was supported by 
GSRP grant NGW-51381. N.F. was partially supported 
by CNPq and CAPES. 



REFERENCES 



Bennett, C. L. 1996, ApJ, 464, LI 

Bond, J. R. in Cosmology and Large Scale Structure, Les Houches 
Section LX, ed. R. Schaeffer, J. Silk, M. Spiro, and J. Zinn-Justin 
469 (Elsevier, Amsterdam, 1996), p469. 

Bond, J. R., Efstathiou, C, & Tegmark, M. 1997, MNRAS, 291, 
L33 

Bunn, E. F. et al. 1994, ApJ, 432, 1.75 
Coble, K. A. 1999, |astro-ph/991141£ 

de Oliveira-Costa, A., Uevlin, IVL, Herbig, T., Miller, A. D., Netterfield, 
C. B., Page, L. A., & Tegmark, M. 1998, ApJL, 509, L77 

Devlin, M., de Oliveira-Costa, A., Herbig, T., Miller, A. D., Netterfield, 
C. B., Page, L. A., & Tegmark, M. 1998, ApJL, 509, L69 

Fisher, K. B. et al. 1995, MNRAS, 272, 885 

Herbig, T., de Oliveira-Costa, A., Devlin, M., Miller, A. D., Page, 
L. A., & Tegmark, M. 1998, ApJL, 509, L73 






Fig. 6. — The top row shows a S/N eigenmodes number 1, 3, 10, 100, 300 and 3082 for the 43 GHz 7UMi map, which are seen to probe 
successively smaller (and noisier) scales. The second row shows cumulative sums up to these mode numbers, with each mode weighted by its 
expansion coefficient so that the total sum recovers the original map (right). The last row shows the analogous running sums with the i*'' 
coefficient rescaled by + 1), so that the total sum recovers the Wiener- filtered map (right). 



, N. & Sillf, J 1997, Nnhm, 386, 37 
preprint pstro-ph/9611212| . 
mionkowsKi, ivi., ivosowsKv, 



ivosowsKy, A & Spergel, D. N. 



Hu, W., Sugiyama, N. & Sill 
Jones, M. E. 1997, 
Jungman, G.. Kamionkows' 

1996, Phys. Rev. D, 54, 1332 
Lahav, O. et al. 1994, ApJ, 423, L93 
Netterfield, C. B. et al. 1997, ApJ, 474, 47 
Rybicki, G. B. & Pvooj. \Ar tt i qq9 J 398, 169 
Staren, J. et al. 2000, Ustro-ph/9912212 ("S2000") 
Tegmark, M. 1997a, ApJL, 48U, Lb I ('*r97a") 
Tegmark, M. 1997b, Phys. Rev. D, 55, 5895 
Tegmark, M. 1997c, Phys. Rev. D, 56, 4514 ("T97c") 
Tegmark, M., de Oliveira-Costa, A., Devlin, M. J., Netterfield, 

B, Page, L. & Wollack, E. J. 1997, ApJL, 474, L77 ("T97d") 
White, M. & Bunn, E.F. 1995, ApJ, 443, L53 

Wiener, N. 1949, Extrapolation and Smoothing of Stationary Time 
Series (NY: Wiley). 



C. 



Wright, E. L. 1996, preprint astro-ph/9612006 
Zaldarriaga, M., Spergel, D.,' fc SfelJ^k, U. 199 ' < , ApJ, 488, 1 
Zaroubi, S. et al. 1995, ApJ, 449, 446 



